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Abstract We have studied the relationship between the solar-wind speed [V] 
and the coronal magnetic-field properties (a fiux expansion factor [/] and photo- 
spheric magnetic-field strength [Bg]) at all latitudes using data of interplanetary 
scintillation and solar magnetic field obtained for 24 years from 1986 to 2009. 
Using a cross-correlation analyses, we verified that V is inversely proportional to 
/ and found that V tends to increase with Bg if / is the same. As a consequence, 
we find that V has extremely good linear correlation with Bg/f. However, 
this linear relation of V and Bs/f cannot be used for predicting the solar- 
wind velocity without information on the solar-wind mass flux. We discuss why 
the inverse relation between V and / has been successfully used for solar-wind 
velocity prediction, even though it does not explicitly include the mass flux and 
magnetic-field strength, which are important physical parameters for solar-wind 
acceleration. 

Keywords: Solar Wind; Radio Scintillation; Corona; Coronal Holes; Magnetic 
fields. Photosphere 


1. Introduction 

Wang and Sheele;^ (|199()1 I1991|l showed that the solar-wind velocity [V] cor¬ 
related inversely with the expansion factor [/] of a magnetic-flux tube in the 
corona, and this inverse relation has been improved by |Arge and Pizzo| (|2()()()H 
for solar-wind prediction as U = 267.5-1-410//^^® (hereafter the WSA equation). 
|Wang, Hawley, and Sheeley] (|1996l) explained this inverse relation with observa¬ 
tional data by the difference of /-dependences of the total-energy flux density at 
the sun and the ion fiux density at the sun: although both parameters increase 
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with /, the increase of the flux density is steeper than that of the total energy 
flux. 

|Fisk, Schwadron, and Zurbuchen] (|1999l) and IFlskl (|2003l) proposed the solar- 
wind acceleration model (hereafter the FSK equation), 


GMq 

2 piUi R© 


( 1 ) 


where V{, Pi, piUi, G, M©, and R© are the final solar-wind speed (e.g. at 1 
AU), the Poyntlng flux at coronal base, solar-wind mass flux at coronal base, 
the gravitational constant, the solar mass, and the solar radius, respectively. 
The FSK equation Indicates the final solar-wind velocity Is determined by two 
parameters, mass flux released from the loops by reconnection with open field 
lines and the Poynting flux at coronal base. 

ISuzukil ((20021 12006(1 introduced the dissipation mechanism of MHD waves 
through shock process to the additional acceleration process of the solar wind, 
and proposed the following equation (hereafter the SZK equation). 


R| Bs 
^^ipvr‘^)iAv f 


{SBa6va)s + ^RTc - , (2) 

7-1 R© J 



where R, 7 , and Bs are, respectively, the gas constant, ratio of specific heats, and 
photospherlc magnetic-field strength. —{6B±6v±)s Indicates the Poyntlng flux of 
the sheared magnetic field, corresponding to Alfveh wave {—{SB±Sv±)s > 0 for 
an outgoing component). Therefore, the SZK equation is a hybrid of WSA and 
FSK models. These three equations have the following similarities: 1) in WSA 
and SZK equations, V is proportional to 1//, 11) In the WSA and SZK equations, 
the /-dependent term works as an additional acceleration, iil) in the FSK and 
SZK equations, the Poynting flux Is used as acceleration energy. Iv) Another 
difference Is the i?s-dependence of V ; although the SZK equation includes Bs 
explicitly, the WSA and FSK equations are Independent of Bs- 

The WSA equation has been widely and successfully used for solar-wind 
velocity prediction because of its simplicity. However, the physical parameters 
necessary for the solar-wind acceleration process, such as mass flux and energy 
are dropped. In order to elucidate the solar-wind acceleration process, we must 
make a model which includes all necessary corona physical parameters explicitly, 
such as the FSK and SZK equations. From this viewpoint, Bs will be one of 
the key parameters. [Wang[ ((2010(1 explained why the WSA equation can drop 
Bs and the mass flux: both the mass flux and energy flux at the Sun tend to 
increase approximately linearly with increasing Bs, but Bs has no systematic 
relation with V at the Earth: accordingly, the Independence of V on Bs causes 
no inclusion of the mass flux and the energy flux In the equation. This means 
that If we found a Bs-V relation, then, necessarily i?s-related parameters of 
the mass flux and the energy flux should be Included. First of all, we will verify 
the WSA equation in Section [31 In Sections 0] and |S] we will examine the Bs-V 
relation and demonstrate that the SZK equation Is the model preferable over the 
others from the viewpoint of solar-wind physics. In Section |6l however, we will 
show that the WSA equation works better than the SZK equation for solar-wind 
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prediction, and explain the reason why Bg and the mass flux can be eliminated 
from the SZK equation. 

This study was began bv IHirano et al.l (|2()()fll) with data obtained from 1995 
to 1996. Analyses were made of large-scale coronal holes (CHs), which consist 
of a polar CH, a distinct equatorial- and mid-latitude CH, and a low-latitude 
polar CH extension. It was revealed that V and Bs/ f have a remarkably high 
correlation, with a correlation coefficient of 0.88 ( |Kojima et al.\ 120041 I2007p . 
Based on their study, here we extend their analyses of data obtained for 24 years 
from 1986 to 2009, including two solar minimum phases. In this study, we were 
particularly interested in comparing the two solar minimum phases of 1994-1996 
and 2007 - 2009, because the latter minimum phase showed peculiar aspects com¬ 
pared with the previous one le.g.. lMcComas et aL[l2008llTokumaru et aLll2009l 
|von Steiger and Zurbuchen| 120111 |Wang and Sheele^ I2013p . Different kinds of 
CHs, such as a polar CH, a middle latitude CH, and a CH associated with an 
active region, have different magnetic-field properties. However, most CHs are at 
middle and high latitudes where spacecraft cannot approach, with the exception 
of Ulysses. Therefore, in this study, we used tomographic interplanetary scintil¬ 
lation (IPS) measurements, which can derive a solar-wind velocity map over all 
latitude ranges ( |Kojima et al.[[r998p . The data and the analysis method will be 
explained in Section [21 


2. Data and Analysis 


In this study we used three data sets of H, Bs and /. We used the solar-wind 
velocities measured at Solar-Terrestrial Environment Laboratory (STEL) using 
the IPS method because the IPS measurements can observe the solar wind at all 
latitude ranges in a relatively short time, and the observations cover several solar 
cycles ( jKojima and Kakinuma[ri990p . However, the solar-wind velocities that are 
estimated directly from the IPS measurement involve the line-of-sight integration 
effect. Therefore, we used the tomographic analysis technique (Kojima et al. 


I1998p . which can deconvolve the line-of-sight integration effects and retrieve the 
intrinsic solar-wind velocity. Using this technique, we derived velocity distribu¬ 
tion maps of latitude and longitude (V-map) on the source surface at 2.5 R 0 
(Figure [T|). Since IPS tomographic analysis requires good coverage of line-of- 
sight projection on the source surface, we used IPS data obtained during three 
Carrington rotations to derive one V-map. Therefore, each V-map represents 
the average structure during three Carrington rotations. This data-window of 
three rotations was shifted by one rotation to derive the map of the next rotation. 
Thus, we obtained about ten maps each year except for the snowy winter seasons 
when the IPS observations at STEL were interrupted. The analyses were made 
for the years from 1986 to 2009 over two solar cycles. 

In order to obtain magnetic-field properties of Hs and /, we used the solar 
magnetic field data observed at the National Solar Observatory/Kitt Peak, and 
calculated the coronal magnetic field configuration by the potential-field-source- 
surface (PFSS) model ( jSchatten, Wilcox, and NessllI969irAltschuler and Newkirkl j 
II969I) . The PFSS analysis code was developed by Hakamada and Kojimaj (112211 
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Magnetic Connection between V and Bs 



Figure 1. A V-map on the source surface at 2.5 R© is derived by the IPS-tomography 
method. Potential magnetic-field lines connect the V-map and the photosphere to obtain the 
flux expansion factor and Bs . 


and was used for mapping a solar-wind structure obtained by IPS to its counter¬ 
part on the photosphere (jOhmi et fll.ll2()()11|Kojima et al. I2()()4|) . We traced back 
the velocity map from the source surface to the photosphere along the potential- 
field line to find the source region of the flow on the photosphere and obtain the 
photospheric magnetic-field strength (Figure[T]). The flux expansion factor [/] 
is defined by the ratio between the magnetic-field strength on the photosphere 
and that on the source surface normalized by (i?ss/R-0)^ : where i?ss is the 
radius of the source surface. The flux expansion factors and the photospheric 
magnetic-field strength were averaged from five solar-wind velocity ranges of 
300-400, 400-500, 500-600, 600-700, and >700 kms“^. Thus, the averages 
of f{Vi), and Bs{Vi) obtained for each Carrington rotation were averaged again 
over a year, from which we derived the year-averaged parameters /, Bs-, and the 
ratio between these two parameters [Bs/f] as functions of V. 


3. Flux Expansion Factor [/] 

Figure [2] shows the dependence of velocity on the flux-expansion factor for the 
years 1986, 1987, 1989, 1993, 1995, 1997, 2000, 2005, and 2009. Triangles with 
solid error bars represent averages of all data in each year. We chose as examples 
the data of nine years among the 24 years, including the solar minimum and solar 
maximum phases. Solid curved lines represent the best-fit solution for 


V = Vo+a/f, (3) 

while the dashed lines represent the best-fit solution averaged for all data during 
24 years given by 

V = 250 + 932//°'30 (4) 
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Figure 2. Dependence of velocity on the flux-expansion factor. Triangles are averages of all 
data in each year. Solid circles are data with weak magnetic-field strength (<5 [G]) and small 
expansion factors (<100), while open circles are those with strong magnetic-field strength 
(>50 [G]) and a large expansion factor (>1000). A solid curved line is the best-fit solution for 
Equation 0. while a dashed line is the average best-fit solution for the 24 years of data from 
year 1986 to 2009. 


The inverse relation between velocity and the flux-expansion factor was origi- 
nally found by Wang and Sheele3^(|199()lll99ip . and its inverse equation has been 
improved by |Arge and Pizz^(l2()()()|) for solar-wind prediction. The original Wang 
and Sheeley equation was based upon the empirical correlation between the ve¬ 
locities observed at the Earth and the flux-expansion factor derived with the data 
measured at the Wilcox Solar Observatory. On the other hand, |Arge and Pizzo| 
(|2000l) propagated the solar-wind predicted on the source surface using Equation 
([3|) to the Earth, and then made iterations changing parameters in the equation 
until a reasonably good agreement was obtained between in-situ observations 
and predicted velocities at the Earth. Although Equation (|1|) can be applied for 
all latitudes, the WSA equation should be applied for the observations near the 
solar Equator because it was derived by making an iterative comparison with 
the low-latitude in-situ data. The numerators of the second term in Equation (j3]) 
and the WSA equation exibit a two-fold difference. This was caused mainly by 
the difference in the magnetic-field data, and therefore Equation 0 should be 
applied to Kitt Peak data, which were measured with a higher resolution than 
those of the Wilcox Solar Observatory. 

Figure [3] shows yearly variations in the expansion factors derived for the five 
solar-wind velocity ranges of 300-400, 400-500, 500-600, 600-700, and > 
700 kms“^. We compared the flux expansions between the last solar minimum 
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YEAR 

Figure 3. Yearly variation in the expansion factors for the five solar-wind velocity ranges of 
300-400, 400-500, 500-600, 600-700, and > 700 kms“^, respectively. 


Table 1. Average expansion factors for the five solar-wind velocity ranges of 
300-400, 400-500, 500-600, 600-700, and > 700 kms“^, respectively. The 
second row is the average through all years [/aii], and the third and fourth rows 
are averages for three-year periods: 1994-1996 [/i] and 2007-2009 [/ 2 ]. The last 
row shows the ratio between two periods. 


V [kms 1] 

300-400 

400-500 

500-600 

600-700 

> 700 

/all 

373±242 

250±151 

31±28 

11±3 

6±2 

h 

221±159 

195±79 

24±2 

12±1 

7±0.4 

/2 

276±21 

167±63 

35±11 

14±0.4 

10±0.5 

/ 2//1 

1.25±0.58 

0.86±0.65 

1.46±0.22 

1.17±0.08 

1.43±0.05 


(2007-2009) and the previous minimum (1994-1996) in Table[TJ The differences 
were less than one-cr variances of the all-year averages, but between two minimum 
phases, there were small but statistically significant differences in medium- and 
high-velocity winds. In the last solar minimum, / increased 1.1 to 1.4 times, 
and this increase was similar to that reported by |Wang and Sheel^ H2013p . 
who compared Ulysses polar passes in 1993-1996 and 2006-2008, finding an 
increase of 1.2 times and a decrease in the median wind velocity by 3% (from 
740 kms“^). 

In order to compare this V-f correlation with two other correlations of V- 
Bs and V-Bs/f in the following sections, we selected two kinds of solar wind. 
One originates from the region with weak magnetic-field strength (<5 [G]) and 
a small expansion factor (<100) (solid circles with dotted error bars in Figure 
[2j hereafter referred to SWl), and the other one from the region with strong 
magnetic-field strength (>50 [G]) and a large expansion factor (>1000) (open 
circles with dotted error bars, hereafter referred to as SW2). We note in Figure 
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Figure 4. Dependence of velocity on the photospheric magnetic-field strength. Triangles are 
averages of all data in each year. Solid circles are data with weak magnetic-field strength (<5 
[G]) and a small expansion factor (<100), while open circles are data with strong magnetic-field 
strength (>50 [G]) and a large expansion factor (>1000). 


[2] that these data not only largely deviate from the average group, but are also 
less dependent on /. 


4. Photospheric Magnetic Field Strength [Ss] 

Figure H] shows examples of the dependence of velocity on the photospheric 
magnetic-held strength for the years 1986, 1987, 1989, 1993, 1995, 1997, 2000, 
2005, and 2009. Triangles with solid error bars are averages of all data in each 
year. An inverse relation between V and i?s exists, but not obviously in the 
three solar minima (1986, 1995, and 2009). Figure [5] shows yearly variation in 
the magnetic-held strength derived for the hve solar-wind velocity ranges of 
300-400, 400-500, 500-600, 600-700, and > 700 kms“^. The average Bq 
during the three years around the two solar minimum phases of 1994 to 1996 
and 2007 to 2009 are shown in Table [21 hr the solar-minimum phase, Bs has a 
very small range of 10.3-12.5 [G] and 6.9-8.2 [G]. This is one of the reasons 
why the inverse relations between V and Bs (triangle) are not obvious in the 
three solar minima. Comparing the two minimum phases, the held strength of 
the last minimum phase decreased about 30 % (Table |2|) of that of the previous 
minimum. [Wang and Sheele^ (|2013l) also reported a 31 % decrease from 10.6 to 
7.3 [G] in footpoint held strength of the polar wind. 
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Figure 5. Yearly variation of the magnetic-field strength of photospheric source regions for 
the five solar-wind velocity ranges same as Figure [3| 


Table 2. Average photospheric magnetic-field strength for the five so¬ 
lar-wind velocity ranges same as Table ^ during three-year periods: 
1994-1996 [Bi] and 2007 — 2009 [B 2 ]. The last row shows the ratio of the 
average magnetic-field strength between two periods. 


V [kms 

300-400 

400-500 

500-600 

600-700 

> 700 

Bi [G] 

B2 [G] 
B 2 IB 1 

12.5±1.7 

6.9±0.5 

0.6±0.3 

11.6±3.5 

7.6±0.6 

0.7±0.5 

10.2±1.0 

7.3±0.6 

0.7±0.2 

10.3±0.2 

7.5±0.4 

0.7±0.1 

10.3±0.4 

8.2±0.6 

O.SiO.l 


The FSK equation indicates that the solar-wind velocity at the Earth should 
have positive dependence on the photospheric-field strength. However, we can- 
not see this tendency in Figure SI and furthermore |Wang, Hawley, and Sheeley] 
(|1996|1 . Wang, Robbrecht, and Sheeley (|2009|l . and | Wang and Sheeley] ((20T^ re¬ 
ported that there appears to be almost no systematic relation between them. 
However, the SZK equation indicates that the additional acceleration term has 
a positive dependence on Bs. Since this term includes not only Bs but also 
some other physical parameters such as / and the mass flux, we examined 
the dependence on Bs, classifying the data by the expansion factor [/]. Data 
are classified into three expansion factors of 5-11, 11-25, and 25-78. Figure 
SI shows as examples three years: 1986, 1989, and 1992, which include solar 
minimum and maximum phases. The figures show that velocity increases with 
Bs when the flux expansion is small, and the dependency becomes less obvious 
as the flux expansion increases. This behavior can be understood by the first 
term in Equation ©; as / becomes large, the first term, which includes Bs/f, 
becomes less effective for determining the final velocity. 
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Figure 6. Dependence of velocity on the photospheric magnetic-field strength. Data are 
classified into three expansion-factor bins, 5 — 11, 11-25, and 25-78. The three years of 1986, 
1989, and 1992 are illustrated. 


In FigurelHwe showed two special kinds of solar wind (solid and open circles). 
Solid circles are the SWl with weak magnetic-field strength (<5 [G]) and a small 
expansion factor (<100), and open circles are SW2 with strong magnetic-field 
strength (>50 [G]) and a large expansion factor (>1000). We note that SWl 
(solid circle) have a positive relation between V and Bs, unlike the data of the 
total average (triangles). 


5. Ratio of Magnetic Field Strength and Flux Expansion Factor 

IBs/f] 

Since the velocity is dependent not only on the inverse of the flux expansion 
factor [1//] but also on the photospheric magnetic-field strength [Bs], we here 
combined these two parameters as Bs //to examine the relation with velocity. 
Practically the parameter Bs/f can be replaced by the magnetic field at the 
source surface [Bss] in the PFSS analysis. Figure [7] shows as examples the 
dependence of velocity on Bs/f for the years 1986, 1987, 1989, 1993, 1995, 
1997, 2000, 2005, and 2009. According to the SZK equation, the correlations are 


SOLA: Fujiki_etal2014SolPhys_color.tex; 15 July 2015; 0:17; p. 9 




















K. Fujiki et ai 


1986 



1987 

- 


' 

£> 

/ / 



. . . / ./ 


.Q 

r/ 

- 


■A/ ■ 


O 

1 / 


O- 

r# • 

■ 

(■>i^ 



/ 


, 


1995 

- 


' 


/ X 



■ ■ ■ 


• 


■ 

- 

X 

- 


■ ■ 




, 


2000 


2005 


E 



1989 



2009 



Figure 7. Dependence of velocity on the ratio of photospheric magnetic-field strength and a 
fiux-expansion factor. Triangles are averages of all data in a year. Solid circles are SWl with 
weak magnetic-field strength (<5 [G]) and a small expansion factor (<100), and open circles 
are SW2 with strong magnetic-field strength (>50 [G]) and a large expansion factor (>1000). 
Solid lines are the best-fit regression lines, while dashed lines are the average for the 24 years 
of data from 1986 to 2009. 


shown for V^. Solid lines represent the best-fit line for the equation of 

V^ = Vo^ + a^, (5) 

while dashed ones represent 

^2 = 2472 -^ 613 ^^, ( 6 ) 

which was obtained by averaging 24 years of data. Although the best-fit line 
changes year by year, had a fairly good linear dependence on Bg/f through¬ 
out the 24 years. 

In order to see the features of the correlation V-Bs/f compared with the 
other two correlations of V-f and V-Bs, we note the two special kinds of solar 
wind of SWl and SW2 shown by open and closed circles in Figure!?] Although 
these kinds of solar wind largely deviated from the major average groups in / 
(Figure [2) and Bs (Figure Sj) correlation diagrams, the SWl data fit well on the 
linear relation with Bs/f. However, SW2 does not have an obvious dependence 
on Bs/f, having a small value of Bs/f. 
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YFAR 

Figure 8. Yearly variation of Bq/ f for the five solar-wind velocity ranges same as Figure 0 


Table 3. Average Bq/ f for the five solar-wind velocity ranges same as Table ^during three-yearl 
periods: 1994-1996 [{Bq/ f)i] and 2007-2009 [{Bq/ f)2]- The last row shows the ratio between two 
periods. | 


V [km s ^ ] 

300-400 

400-500 

500-600 

600-700 

> 700 

{Bs/f)i 

{Bs/f )2 

{Bs/f) 2 /{Bs/f)l 

0.176±0.072 

0.131±0.027 

0.74±0.62 

0.297±0.056 

0.187±0.036 

0.63±0.43 

0.599±0.142 

0.383±0.056 

0.64±0.44 

1.101±0.087 

0.620±0.055 

0.56±0.21 

1.778±0.019 

0.906±0.058 

0.51±0.13 


Figure [HI shows the yearly variation of Bs/f derived for the five solar-wind 
velocity ranges of 300-400, 400-500, 500-600, 600-700, and > 700 kms“^, 
and Table [3] compares variations between the two solar minimum phases. Bg/f 
in the last minimum phase decreased from that in the previous phase, and high- 
velocity winds (600-700 and > 700 kms“^) especially decreased by as much as 
50 %. The major percent of this decrease was caused by the 30 % decrease of Bs 
(Tabled]) and the 10 to 40% increase of / (Table [T|). 


6. Discussion 

We investigated the relation between the solar-wind velocity and the coronal 
magnetic parameters from 1986 to 2009. From the cross-correlation analysis, 
we verified that V has an inverse relation with / (Equation ([1|)) as reported 
by jWang and Sheeleyj (|1990111991|) . There are reports that both the mass and 
energy-flux densities tend to increase approximately linearly with increasing 
magnetic-field strength, but there appears to be almost no systematic relation 
between photospheric magnetic-field strength and the solar-wind velocity at the 
Earth ( jWang, Hawley, and Sheel^ 119961 [Wang, Robbrecht, and Sheeley|l200T?l 
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1990 1995 2000 2005 2010 


YEAR 

Figure 9. Differences in the best-fit solutions between the solid line and the dotted line in 
Figure [2] The discrepancies are shown for four different expansion factors of 2663, 489, 160, 
and 69, which give velocities of 400, 500, 600, and 700 kms“^from Equation 0. respectively. 


|Wang and Sheele^ I2013|l . However, we found a different relation between pho- 
tospheric magnetic-field strength and the solar-wind velocity. Investigating the 
V-Bs relations for the solar wind binned by the flux expansion factor, we found 
that the higher velocity winds flow out of stronger field regions in each bin. 
Combining these two relations of V-f and V-Bg , we revealed that the Bg/f has 
an extremely good linear correlation with the square of the velocity (Equation 
El- This linear correlation also holds for the solar wind from CHs that have 
weak magnetic held and a small expansion factor, although they substantially 
deviated from the V~f correlation. 

We compared Equations and ([H]), which are better models for making 
solar-wind predictions through the whole solar cycle. Since these equations were 
derived from averaging 24 years of data, we examined discrepancies in them 
from the equations obtained for each year (discrepancies between the dashed 
line and the solid line in Figures [5] and [7|). Figure O examines Equation (|3|). 
The discrepancies are shown for four different expansion factors of 426, 78, 25, 
and 11 which give velocities of 400, 500, 600, and 700 kms“^from Equation 
O, respectively. The average discrepancies over the 24 years are —11 ± 23, 
—22 ± 40, —30 ± 48, and —30 ± 48 kms“^at each of the four expansion factors. 
This means that a prediction over the whole solar cycle can be made accurately 
using the single equation. In particular, the prediction of low- and medium- 
velocity solar winds, which are usually observed at the solar Equator, can be 
made with accuracy better than 50 kms“^. 

Figure [TUI examines Equation ([51). Discrepancies between dashed lines and 
solid lines in Figure[7]are shown for three different values of Bg/f: 0.264, 0.503, 
and 0.796, which give velocities of 400, 500, and 600 kms“^ from Equation 
®. Figure [10] indicates that this single equation cannot be applied for solar- 
wind prediction through the whole solar cycle. This is because not only Bg/f 
but also the parameter a change significantly from year to year. This result 
reflects following two facts: i) the parameter Bg/f which is exactly same as the 
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Figure 10. Difference in the best-fit solutions between the solid line and the dotted line in 
Figure [71 The discrepancies are shown for three values of B^/f of 0.068, 0.102, and 0.144, 
which give velocities of 400, 500, and 600 from Equation respectively. 


V^=V(,2 + a*(B/f) 



magnetic field strength on the source surface in this analysis determines the 
global magnetic structure in interplanetary space and depends mainly on solar 
magnetic dipole moment which changes significantly through solar cycle, and ii) 
the range of solar wind velocity hardly changes through solar cycle. 

It has been reported that the polar solar wind decreased in velocity in 2006 - 
2008 compared with that observed in 1993-1996, but the decrease was as small 
as 3% (from 763 to 740 kms“^) (IMcComas et al.[ 120081 |Wang and Sheele^ 
120131) . On the other hand, the changes of Bs, / and Bg/f between the two 
minimum phases were not small. In particular, the change of the parameter a in 
Equation ([5|) was as large as 200% between the two minimum phases ('Figure ITT] 
and Table 0]). The less significant change of the polar high-velocity wind between 
two minimum phases can be understood by the second term of Equation 
aBs/f, which gives additional acceleration to the thermal-driven solar wind. 
Since a increased 200% and Bs/f decreased about 50% (Tabled, the change 
of the additional acceleration is small. 
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Table 4. Comparison of averaged parameters Vq and a in Equa¬ 
tion © between two solar minimum phases, The third column 
represents avreaged B^/f for V = 750 kms~^. The last row 
shows the ratio between the two periods. 


period 

Vo 

a 

Bs/f [750 kms-i] 

1994-1996 

321±20 

535^ ±57^ 

1.778±0.019 

2007-2009 

237±43 

7822±2552 

0.906±0.058 

2007-2009 

1994-1996 

0.74±0.26 

2.13±0.05 

0.51±0.13 


Table 5. Average solar-wind density at 1 AU for the five solar-wind velocity 
ranges same as Table[T]during three-year periods: 1994-1996 [A^i] and 2007-2009 
[N 2 ] - The last row shows the ratio between the two periods. 


V [kms 

300-400 

400-500 

500-600 

600-700 

> 700 

Ni [cm“®] 
N 2 [cm“^] 
N 2 /N 1 

11.13±6.29 

6.77±4.84 

0.61±3.85 

6.70±4.30 

4.11±3.01 

0.61±2.68 

4.29±2.17 

2.96±1.74 

0.69±1.55 

3.49±1.34 

2.58±0.89 

0.74±1.02 

3.19±2.32 

2.01±0.61 

0.63±1.47 


The square of the first term in Equation (0, Vq , corresponds to 2[7/(7 — 
l)Rrc —(GM 0 /R 0 )] in the SZK equation. Therefore, tQ^/2+GM0/i?0 is propor¬ 
tional to the coronal temperature [Tc], and the Vq decrease of 25 % between two 
minimum phases (Table0|) corresponds to a 6 % decrease of coronal temperature. 
This is the same freeze-in temperature decrease derived from the ion charge ratio 
of solar wind by |von Steiger and Zurbuchen] (|20lT]) . 

According to the SZK equation, a varies in proportion to the inverse of mass 
flux. Although the solar-wind velocity did not change signihcantly between the 
two minimum phases, solar-wind densities provided by the NASA/GSFC’s Space 
Physics Data Facility’s OMNIweb service decreased from 30 to 40% (Table 
[5]). IMcComas et al.\ (I2008D also reported that the fast solar wind from both 
large polar CHs was significantly less dense (17%) than during the previous 
solar minimum. Therefore, half of the a increase could be attributed to this 
density decrease. Since a also consists of the Poynting flux from the photo¬ 
sphere another half of the a increase may have been caused from 

an increase of the Poynting flux. From this viewpoint, it is worth noting that 
|Wang and Sheel^ (I199ip commented that if Alfven waves were responsible for 
accelerating the solar wind, and the wave energy flux is roughly constant, the 
mean-square wave velocity [((^Ug)] and the field strength [Bg] at the coronal base 
may be inversely correlated. 

Making cross-correlation analyses between the solar-wind data measured by 
Advanced Composition Explorer (ACE) spacecraft and the photospheric magnetic-| 
held strength, Wang and Sheeley (2013) obtained the following equation between 
the mass flux density at the coronal base and the footpoint held strength: 
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(7) 


where i?s is a footpoint field strength. Therefore, if we compare the solar-wind 
streams which have the same flux expansion factor, Bs/ueVe is a constant value 
that does not depend on solar activity. In fact, Tables [2] and 0 show that the 
changes of Bg and ngVe between the two minimum phases are the same amount 
of 30%. Consequently, the mass flux and the magnetic-field strength disappear 
from the additional acceleration term in the SZK equation, becoming const, x 
{5B±Sv±)s/f, which is similar to the WSA equation except for the Poynting 
flux. This equation can be made equivalent to the WSA equation, by making the 
Poynting flux less dependent on the solar activity. However, we discussed above 
that the Poynting vector has to be about 1.4 times larger in the last minimum 
phase to make the additional acceleration term in the SZK equation less sensitive 
to solar activity. This contradiction must be investigated in a further study. 


7. Conclusions 

We investigated experimentally the relation between the solar-wind velocity and 
the coronal magnetic parameters by comparing the WSA and SZK equations. We 
verified that V is inversely proportional to / and found that V tends to increase 
with Bs if / is the same. Consequently, we revealed that V has an extremely 
good linear correlation with Bg/f as indicated by the SZK equation. However, we 
should note that the SZK equation {V-Bs/f relation) cannot be used for solar- 
wind prediction through the whole solar cycle. The slope of its linear relation 
changes year by year largely because it explicitly includes a mass flux. On the 
other hand, the WSA equation [V-f relation] is a good prediction model, which 
can be used through the whole solar cycle, although it includes neither the mass 
flux nor Bs which play important physical roles for the additional acceleration 
of the solar wind. This is because the mass flux and magnetic-held strength can 
be canceled in the SZK equation using their linear relation. 

The slope of the linear relation, V-Bs/f, shows a solar-cycle variation. Es¬ 
pecially we focused on the comparison between the two solar-minimum phases, 
1994-1996 and 2007-2009. Consequently we found that the slope is two times 
larger in the very quiet solar minimum (2007-2009) in comparison with that in 
1994-1996. This difference in solar-minimum phases is probably explained by a 
long-term variation of the solar-cycle activity. However, the slope also tends to 
increase during solar active phases, which was not discussed in this article; we 
will extend this examination including the whole solar cycle in a further study. 
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